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Abstract 



We have extended the variational perturbative theory based on the back 
ground field method to include the optimized expansion of Okopinska and the 
post Gaussian effective potential of Stansu and Stevenson. This new method 
provides much simpler way to compute the correction terms to the Gausssian 
effective action (or potential). We have also renormalized the effective po- 
tential in 3 + 1 dimensions by introducing appropriate counter terms in the 
lagrangian. 
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I. INTRODUCTION 



The effective action and the effective potential have been used effectively in studying 
various aspects of quantum field theories. Nonperturbative approximation methods for cal- 
culating them have been used to extract essential features of nonlinear quantum field theo- 
ries, such as properties of the ground state and information about possible phase transition 
of the system. 

One of the most convenient nonperturbative approximations to the effective potential is 
the Gaussian variational approximation proposed more than 20 years ago [1]. The corre- 
sponding Gaussian effective potential (GEP) is known to contain one loop, sum of all daisy 
and superdaisy graphs of perturbation theory [2] and leading order in 1/7V expansion. 

Lately some procedures to compute the correction terms to the GEP were proposed [3-6] . 
In our previous paper [6] we have developed a systematic method of the perturbative ex- 
pansion around the Gaussian effective action based on the background field method (BFM). 
Although BFM provides one with the simplest method of computing the perturbative cor- 
rections to the GEP, one can extend the method to give a better approximation. The point 
is that in refs. [5,6] the effective mass (Q) is fixed by the Gaussian gap equation, whereas Q 
is optimized at each stage of approximation in refs. [3,4]. 

The main purpose of the present article is to extend BFM to include various other post 
Gaussian methods developed so far. We shall develop an improved post Gaussian effective 
potential (PGEP) method using BFM, where Q is not fixed preliminary. 

The paper is organized as follows. In Sec. II we briefly present the basic BFM formulas, 
which will be necessary for further developments. In Sec. Ill we derive PGEP using BFM for 
the 4 theory, and we show that with BFM the computation of the higher order corrections is 
very much simplified compared to other methods. In Sec. IV we consider the renormalization 
procedure. The summary and conclusions are given in Sec. V. 

II. BACKGROUND FIELD METHOD FOR SCALAR FIELDS 

For readers convenience we outline here the basic formulas of BFM. It is well known 
that, BFM is the most effective way of doing the so called loop expansion [7] allowing us 
to study one - particle irreducible (1PI) Feynman diagrams. The generating functional for 
disconnected graphs 1 is defined by 

Z[j] = J D0exp{*(S[0] +j<f>} = (0\0y (2.1) 

where j(x) is the external source and S[4>] is the classical action: 

S[<t>] = J d 4 xC[<l>(x),d^(x)}. (2.2) 

In BFM an analogous quantity, in which the classical action is written as a function of the 
field plus an arbitrary background field B, is defined by 



Here and below we use integral convention, e.g. j(f> = / d xj(x)<f>(x) . 
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Z[j,B] = J V<Pex P {i(S{ ( f ) + B] + #}, 



(2.3) 



where Z[j,B] depends both on the conventional source j and on the background field B. 
Also by analogy with the conventional generator of connected graphs, 

W\j] = -i\nZ\j], (2.4) 

we define a new generating functional for connected Green's functions, 

W\j,B] = -i\nZ\j,B]. (2.5) 

In the conventional approach the vacuum expectation value of the field operator in the 
presence of external source is defined by 

Mx) = {<f>(x))' = ^. (2.6) 

In BFM it is defined as 

5j(x) 

and the background field effective action is defined by 

T[^,B]=W[j,B]-j^. (2.8) 

Now it can be easily shown [7] that the conventional effective action can be calculated by 
the following equation: 

r[0o] = f[O,5]. (2.9) 

The background field effective action T[4>q, B] is just a conventional effective action computed 
in the presence of background field B. It therefore consists of all IP I graphs contributing 
to Green's functions. Since 1PI Green's functions are generated by taking derivatives of 
r[0Oj-B] with respect to 0o it would generate 1PI Green's functions in the presence of 
the background field B. Now r[0, B] has no dependence on O) so it generates no graphs 
with external lines. Instead, r[0, .B] is the sum all 1PI vacuum graphs in the presence 
of the B field. The advantage of the background field method is that many diagrams in 
the conventional perturbation method are amalgamated into one diagram in BFM due to a 
new definition of the propagator function [7]. This method can be modified to extract the 
Gaussian part of the effective action out of the perturbative action part [6]. We will now 
show that BFM can be used to compute the perturbative expansion of the effective action 
around the Gaussian effective action wihout fixing the effective mass. 
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III. POST - GAUSSIAN EXPANSION OF THE EFFECTIVE POTENTIAL 



For technical convenience we shall work throughout in the Euclidian formalism and start 
with the Lagrangian of the 4 theory: 

£=l0(-9 2 + m2)0 + A o 4 , (3.1) 

where m Q and A are bare mass and bare coupling constant, respectively. Now following 
refs. [3,4] we rewrite (3.1) in an equal form by introducing a new mass parameter Q : 

C = l -<p(-d 2 + fi 2 )0 + l -{ml- Ql)^ + A o 4 . (3.2) 

The generating functional Z[j,B] is then given by Eq. (2.3): 

Z[j,B} = |p0exp{-i(5[0 + J B]+j0)} = Jv^eM-liJd'xjC^ + Bj+m, (3.3) 

where C(<f> + B) can be separated as follows: 
£(0 + 5) = £ o (0) + £i(0,5), 

C = |0(-9 2 + & )<t>, (3.4) 
B) = v (B) + Vl (B)<j> + v 2 (B)(j) 2 + v 3 (B)(j) 3 + v A {B)^ , 



with 



m 2 B 2 

VQ (B) = ^- + \ Q B\ v 1 (B) = B(ml + A\ B 2 ), 



v 2 {B) = |K - «o) + 6A o5 2 , v 3 {B) = 4A £?, ( 3 - 5 ) 
v 4 (B) = A . 

Now performing explicit Gaussian integration in Eq. (3.3) one obtains 

Z\j,B] = exp{-i| d^xC^ -> 5/5j, B)} J T?0exp{-(i| d A xC Q {<P) - 
= [detV- 1 ]-^ expj-i J tPxdty -> 6/6j,B)}exp{hjVj/2} , 

where 

^ = (-^ 2 + ng)4* ■ (3.7) 

In accordance with our previous work [6], we introduce the so called primed derivatives: 



(3.6) 



5 2 

—)' 




5 2 




SfJ 






5 3 

—)' 


= 4 3) = 


5 3 


3/i "Z^ xx^^x 


SfJ 








5 A 

—)' 

8j? 


- ii 4) = 


5 4 


5 2 

- 6hv xx — + 3h 2 V xx 2 



(3.8) 



where R x = J d 4 yV xy j(y), so that 

A^exp{hjVj/2} = h n R n x exp{hjVj/2} . (3.9) 
To isolate the Gaussian approximation we introduce another Green's function: 



exp{-±(v 2 5 2 /5j 2 + v 3 5 3 /5f + v 4 6 4 /5j*)}exp{hjVj/2} 

= N ex V {\a{B)5/5 3 }ex V {-\{v 2 M 2 ) + v 3 A^ + v 4 A^)}ex P {HjGj /2} , 



(3.10) 



where N , a(B) and G are to be determined. By using the definition (3.8) one finds the 
following solution to the Eq. (3.10): 

a(B) = -3v 3 (B)hG xx , 

N (B) = (detG- 1 )- 1 / 2 (dctV- l ) 1 / 2 exp{-v 2 G xx + 3v 4 hG 2 xx }, 1 ' j 

where G xy satisfies the equation: 

G x i = Ky + l2hv 4 G xx S xy . (3.12) 
The last equation looks rather simple in momentum space: 

n 2 = n 2 - i2h\ i (n), (3.13) 

with 

_ f d 4 pexp[-ip(x - y)\ 



r dp ex 
" J (2tt 



(2nWp 2 + n 2 ) ' 



G, y = i d4p : x ±: p{x - y) \ (3.H) 



Jo(ft) 



(2tt) V + Q 2 ) ' 
d A p 



(2tt) 4 (p 2 + Q 2 ) ' 

As to the determinant factors in Eq. (3.11), the first one can be written in a usual form: 

(detG~r i/2 = eM-h(m, /i(n) = | / d4pl ^" 2) . (3- 15 ) 

and the second one will be canceled by the term (det D -1 ) -1 / 2 in Eq. (3.6). Thus substituting 
Eqs. (3.10), (3.11), (3.14) into (3.6) we come to the final expression for Z[j, B] and W\j, B\. 

Z[j, B] = exp{-i[n/ 1 (fi) + v + hv 2 G xx - 3v 4 h 2 G 2 J}I B , 

W[j, B]/h = In % 5] = — - Vft - «2G M + 3v 4 hG 2 xx + In J B , 

with 

7 B = exp{-i( Vl + 3u 3 fcG xx )<W} exp{-i[u 2 i< 2 > + t; 3 i {3) + v 4 i (4) ]}exp {HjGj/2} . (3.17) 
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The effective potential V e // is defined by 



Veft(B) = -^Sl (3.18) 
/ 

when r(_B) is independent of momenta. Using Eqs. (2.5)-(2.9) and (3.16)-(3. 18) one can 
represent V e ff as a sum of the GEP, V G , and a correction AV^: 

V e// (fl) = V G (B) + AVb(S) . (3.19) 

The explicit expression for AV G will be given below. Here we note some points. 

1) Due to the special construction of primed derivatives (3.8), the coefficient of I B in 
Eq. (3.16) gives rise to the Gaussian effective action, so from Eqs. (2.9), (3.16), (3.18) we 
immediately obtain the GEP: 



i 2 G 2 

= hh(Q) + vq(b) + v 2 hi (n) - 3v 4 ti 2 i$(ty 



v G {B) = hh(n) + v (B) + v 2 UG xx - Sv 4 h z Gl x 



hh(n) + \m 2 B 2 + A 5 4 + M^)\\(m 2 - Q 2 ) + 6\ B 2 ] - 3\ h 2 I 2 (Q) 

hh(n) + \m 2 B 2 + a £ 4 + hi (n)[±(m 2 - n 2 ) + 6\ b 2 ] + 3\ h 2 i 2 (n) 



(3.20) 



Here we used the explicit expressions for vq,v 2 ,V4 given by Eqs. (3.5) and (3.13). 

2) The correction to the GEP, that is, AV G) in Eq. (3.19) is given by \nI B - As it was 
explained in [6] the linear term in the exponent in Eq. (3.17) can be omitted, since it does 
not contribute to the effective potential. Thus the corrections to the GEP will be given by 

AV G (B) = -h\nI B = -/Un{exp[-fM( 2 ) + v 3 A^ + v 4 A^)}exp {hjGj /2}}\ J=0 



_ Mn{1 _ 6(v 2 AW + v 3 A® + v A A^) expjhjGj/2} 



=0+ 



5 2 (v 2 A^ + U3 i(3) + V4 AW) 2 exp[hjGj/2] . ( 3 - 21 ) 



+ 2h 2 b=0 



5 3 (v 2 A^ + v 3 A^ + v A A^f exp[hjGj/2] 

67? Ij=0 + • • - } 

= 6AV$\b) + 5 2 AV G 2) {B) + 5 3 AV G 3) (B) + ... . 

Here we have introduced an auxiliary expansion parameter 5(5=1). 

The first order term AVq\b) in this equation will not contribute to the ef- 
fective potential, i.e., AV G \B) = 0, due to the relation, A^exp {hjGj /2}\ j=0 = 
h n R%exp{hjGj/2}\ j=0 = (see Eq. (3.9)). 

3) As to the next term of the Eq. (3.21) of order 5 2 the explicit calculations show that, 
only diagonal terms in the expression (v 2 A^ + v 3 A^ + v A A^) 2 give a nonvanishing result 
at j = 0: 



6 



(A( 2 )) 2 exp{%'Gj/2}|, =0 = 2\h 2 J d*yG 2 xy 

(i( 3 )) 2 exp {hjGj/2}\ j=0 = 3\k 3 J d'yG 3 xy (3.22) 
(iW) 2 exp {hjGj/2}\ j=0 = A\h 4 J d A yG xy . 

Therefore, to 5 2 order we obtain the following expression for the first order correction to the 
Gaussian effective potential, 

av£\b) = -^[v 2 2 i 2 (n) + hv 2 3 i 3 (n) + h 2 v 2 Mm , (3.23) 

where in accordance with ref. [4], the following integrals are introduced: 

J»(n) _ r + ^ fd* P G 2 (p) 



2! - J d y ° xy - J (2tt) 4 ' 
hi®) _ f j4 ^ 3 f d*pd*q G(p)G(q)G(p + q) 



f ,4 r 3 _ [ d-pd-qG{p)G{q)G{p + q) 

J yLr *y ~ J (2vr)8 ' (3.24) 



3! ~J * xy J (2tt) 

_ /• , 4 „ 4 f d A pd*qd A k G{p)G{q)G{k)G{p + q + k] 



d 4 yG 4 2 



4! j ~ * xy j ( 27r p 
G(p) = l/(p 2 + n 2 ). 

4) The expressions like {A^A^ . . .}exp {?j,j'G r j/2}| J= o can be calculated analytically by 
using Mathemathica. The calculations show the following result in 5 3 order 2 : 



A^AfAfexp {hjGj/2}\ j=0 = 8G xy G yz G zx , 
ifij^exp {hjGj/2}\ j=0 = 2AG 2 xy G 2 yz , 
AfAfAfexp {hjGj/2}\ j=0 = 36G xy G 2 yz G zx , 



yz i 

(3.25) 



A x VA y ^A^exp{hjGj/2}\ J=0 = 192G xy G 3 yz G zx , 
A^AfAWexp {%'Gj/2}| J=0 = 216G 2 ^GL , 
A^A y ^A z Vexp{hjGj/2}\ 3=0 = 1728G 2 xy G 2 yz G 2 zx . 

Using these equations in (3.21) we get the second order correction (5 3 terms) to the GEP, 

AV^ (B) = \v\G xy G yz G zx + \2v 2 v A G 2 xy G 2 yz + I8v 2 v 2 G xy G 2 yz G zx + 
+96v 2 v 2 G xy G 3 yz G zx + 108v 2 v 4 G 2 xy G yz G 2 zx + 288v 3 G 2 xy G 2 yz G 2 zx . 

We see that, correction to the GEP in 5 2 , i.e. , AVq\b) consists of only three Feynman 
diagrams shown in Fig. 1 (a), which we call BFM diagrams. In Fig. 1(a) quadratic vertex 
(marked as ^> ) represents 



2 A11 other combinations vanish. 
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v 2 (B) = -(m 2 -n 2 )+6\ B 2 , (3.27) 

On the other hand, AVq\b) includes five diagrams found in ref. [4]. These two additional 
diagrams may be obtained by using (3.13) in (3.27) and further in (3.23). As a result 
one gets all diagrams of ref. [4] shown in Fig. lb. One can see that the diagrams of ref. 
[4] can be obtained from the BFM diagrams by attaching a ring diagram with Io(Q) = 
J <i 4 p/(27r) 4 (p 2 + f2 2 ) to each v 2 vertex of a BFM diagram. 

Six BFM diagrams in 5 3 order given by Eq. (3.26), are presented in Fig. 2a. Now 
using Eqs. (3.13) and (3.27) in (3.26), i.e. , by attaching a ring diagram , denoted by a 
small circle in Fig. 2a to each v 2 vertex, one gets seven additional diagrams as is shown in 
Fig. 2b. After this procedure AVq\b) will be the same as the one given by Okopinska [3], 
represented by 13 diagrams. One may conclude that, the application of BFM technique to 
study corrections to the GEP makes the work much easier than in the formalism used in 
refs. [4,3]. To illustrate this we represent in Fig. 3 nine BFM diagrams including t> 2 vertices 
in 5 4 order. The implementation of a ring diagram, through Eqs. (3.13) and (3.27) gives 
additional 18 diagrams. 

The effective potential in Eq. (3.23) is exactly the same as that given in refs. [5,6], except 
that the parameter fl in Eq. (3.23) is not fixed by the Gaussian part. In fact in accordance 
with the principle of "minimal sensitivity" Q will be fixed by dV e ff/dfl = 0. 

IV. REGULARIZATION AND RENORMALIZATION 

It is well known that, the renormalization of any field theory can be achieved by intro- 
ducing appropriate counter terms to the Lagrangian. In this section we apply this procedure 
to evaluate V e ff- To do this we rewrite the lagrangian (3.1) including counter terms [9,10]: 

C = \(t){-d 2 + m 2 )0 + A0 4 + ]-Bm 2 (t) 2 + CA0 4 - Dm 4 (4.1) 

where m and A represent renormalized mass and renormalized coupling constant, respec- 
tively, and make the following mass substitution: 

m 2 = Q 2 - 5(Q 2 - m 2 ) = n 2 -s x (4.2) 

where 5 is a superficial parameter, well known in the delta expansion method. We then 
rewrite the Lagrangian in the same way as in 5 expansion [9] 

C = C + 5{C - Co) = ^(-d 2 + Q 2 )<p + C mt (4.3) 

where 

Ant = ~¥x<P 2 + + £ cnt 

(4.4) 

£ cnt = \<p 2 B{n 2 - 5 X ) + AC0 4 - D(Q 2 - 5 X ) 2 • 
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The idea that the mass substitution (4.2) should be used not only in the standard mass 
term but also in the counter terms was first recognized in [11] and further developed in [10] 
where the nessesity of the counter term D was also pointed out. 
We assume that the parameters B, C and D are expanded as 

B = B v 5 + B 2 5 + 0(5 3 ) 

C = d5 + C 2 5 2 + 0{5 3 ) (4.5) 
D = D + D 1 5 + D 2 6 2 + O(5 3 ). 

Now using the method outlined in the previous sections one may get the following un- 
renormalized effective potential up to the S 2 order 



where 



\/ e //(0o) = Vb(0o) + A\/ G (0 o ), 

Vg(0 o ) = h(n) - D Q 4 + 5{v + v 2 I (Q) + 3vJ 2 (Q)) 

AVb(0 o ) = -5 2 {\{v 2 + 6vj (n)) 2 i 2 (n) + lv 2 i 3 (n) + \v 2 h{n) 
-(u + u 2 i (n) + 3u 4 i$(n))} , 



(4.6) 



v 2 = -\x + 6A0g + -J— + 6CiA^ , 

v 3 = 4A0o + 4A0 o Ci , v 4 = A + Ad , (4.7) 

(B 2 tt 2 - -Bix)0p 4 nn 4,onn2. n.2 

m = ^ h AC 2 o - -L> 2 S2 + 2L>iii x - M)X , 

« 2 = (£?2 ^ ~^ lX) +6AC 2 ^, M3 = 4A0 o d, « 4 = AC 2 , 

and the divergent integrals J„(f2) calculated in the dimensional regularization method in 
3 + 1 dimensions are brought to the Appendix. 

We shall evaluate the regularized effective potential 3 order by order in powers of 5. 

It is easy to see that, the first divergent term of the GEP , coming from in Eq. 

(4.6) may be compensated by D chosen as: 

32tt 2 e 

As it is seen from Eq. (4.1), possible finite part of the counter term D leads only to shifting 
the effective potential as a whole and may be neglected. Using explicite expressions for 
divergent integrals in Eqs. (4.6) and (4.7)the GEP is given by: 

\/ G (0o) = A(d + 1)0 O 4 + M( 2 )(fi)0o 2 + M^\n) + M(°)(fi) (4.9) 

where 



3 More exactly V eff (<p , fl 2 ) - V ef f (0,m 2 ) 
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n 2 

M(2)ro . 3^A(Q + l)ln(-) [3A(Cl + 1)(7 _ 1)+4gl7r 2 ]n 2 + 4m 2 7r 2 

M {n) = 8^ + 8^ ( 4 - 10 ) 

3^ 2 A(Ci + 1) 

47r 2 e 

and Af(°)(fi) will be given below. From the first term of Eq. (4.9) one may conlude that, 
C\ should be finite. The requrement of finiteness of M^(Q) in Eq. (4.10) leads to : 

fll = iB + «f£22 (4.H) 

2tt 2 € 

where bo, as well as C\ are finite constant to be determined by suitablely chosen renor- 
malization condition. Fortunately, the choice of B\ as in Eq. (4.11) will cancel also the 
logarithmic pole like ln(Q 2 /fi 2 )/e coming from Iq(Q) in Eq. (4.6). The term M^\fl) in Eq. 
(4.9) includes the logarithmic part, 



m 2 

3 A(Gi + lj(irin'( 



A(C 1 + l)(Q 4 m 2 (^)-m 4 ln 2 (^)) 

Xi 

256 vr r 



r(0) (ci\ - 3 



+ T7^{[(3A(C 1 + 1)( 7 - 1) + Ab n 2 - 27r 2 )fi 4 + 4m 2 7r 2 fi 2 ]ln(^) (412) 

2 

- [m 4 (3A(d + 1)( 7 - 1) + 46 tt 2 ) + 2m 4 vr 2 ]ln(^-)} , 

which is now finite, and the nonlogarithmic part 

M(°)(fi) = — ~ m 4 )( 7 - l)[3A(d + 1)( 7 - 1) + 8b 7r 2 ] + 27r 2 (3m 4 + Q 4 ) 
2567T 4 

+ 8m 2 fi 2 7r 2 ( 7 - 1) - 47r 2 7 (m 4 + fi 4 ) ( 413 ) 

, / 4 4n , 16&Q7T 2 , 12A(d + l) n 
+ (m 4 - ir) 2567r 4 L>i H 1 J} . 

e e A 

which becomes finite when D\ is fixed as : 

Dl = -i^-'^+i). (414) 

16 7r 2 £ 64 7r 4 £ 2 

Now, after having removed all the singularities one comes to the following finite effective 
potential in 8 order, 
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'5 2 

V^o) = ^mV+ m + 1)0 O 4 + — >: g7r2 



^ 2 0o 2 pA(d + l)ln(^-) + 3A(d + 1) + 4b 7r 2 



+ 



3A(C 1 + l)[fi 4 ln 2 (^)-m 4 ln 2 (^)] 



256tt 4 



+ T^[(( 3A ( C i + !) + 4b o^ 2 " 2vr 2 )^ 4 + 4m 2 vr 2 fi 2 )ln(^) 



m 2 

m 4 ((3A(d + 1) + 4b 7r 2 ) + 27T 2 )ln( — )] 

H 1 



+ 25^ [( ^ 4 " m4)(7 " 1 )( 3A (° 1 + X) + 8bo7r2) + 27r2 ( 3m4 + Q4 ) 
+8m 2 fi 2 7T 2 ( 7 - 1) - 47r 2 7 (m 4 + Q 4 )] . 



(4.15) 



Note that, V^g(0o) includes two finite constants b and Ci to be determined by the suitably 
chosen renormalization conditions. If we employ the minimal subtraction (MS) scheme 
these constants may be neglected, since in the MS scheme the counter terms are supposed 
to remove only the singularities in e. On the other hand, they may be determined by an 
intermediate renormalization scheme [12]: 

d Xs(M | m 2_ i d %s(<f>°) i A _ Ul6) 

# 2 k=0 " m -°' 4! M U = ~ A -°' (416) 

which is commonly used for fixing the parameters of the effective potential. In this case 
they may naturally depend on ln{m 2 / n 2 ). From the first condition in Eq. (4.16) the finite 
counter term b is easily determined as 

,m 2 



3 A(C 1 + l)(7 + ln(— )-l) 
"4 ^ 



b* = -- A ^-^ , (4.17) 



while C\ can be determined in principle from the second condition of Eq. (4.16) which 
turned out to be a nonlinear equation. We shall come back to this point later. 

Now passing to the regularization of V e g in next order of 5 we rewrite it explicitely as: 

avg(0 o ) = K^mi + [K%\n) + K^tmi + + /r<o)(n) (4.18) 

where 

A[9A(C 1 + l) 2 ( 7 + ln(g)) + 4C 2 , 2 ] g (4ig) 
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and if( 2 )(f2), K(°\£l) will be given below. Finiteness of K^ 4 \il) determines C 2 as 

C 2 = C 2 + 9A( ^ 2 +1)2 (4.20) 
2ir 2 e 



where C 2 is a finite constant. The term K^\Q) is given by: 



32 7T 4 



3 O 2 m 2 

Aln(^){3A(C 1 + l) 2 fi 2 ln(^-) 



32vr 4 > 2 v 1 y V 2 ' (4.21) 

+ [(36A - 27A 7 )Ci 2 + (72A + 4tt 2 - 54A 7 )d + 4tt 2 + 36A - AC 2 tt 2 - 27A 7 ]fi 2 

A 2 fi 2 (d + l) 2 lnfi 
- 4tt 2 (C 1 + l)m 2 } - | — 

Here we see that there is a dangerous pole term propotional to ln(0 2 //x 2 )/e which can be 
removed only by the the condition 

Ci = -1 , (4.22) 

since there is no infinite counter term in (4.21) to compensate this pole. Now the nonloga- 
rithmic term in (4.18), K^ 2 \Q), is simply given by: 

K(2) = I [4g 2 7r 2 + 3AC' 2 (-1+ 7 )P 2 _ 3 AC^ (423) 

8 7T 2 4 7T 2 £ 

and may be regularized by 

B 2 = 6 2 + H (4.24) 

where 6 2 is a finite constant . One can see that the choice Ci = — 1 will remove logarithmic 
pole terms in (Q) also : 

m 2 O 2 



AC 2 [m 4 ln 2 ( — )-n 4 \n\— 
256 7^ 



< j (^) 

Z,dU /I " 



x (4fi 4 6 2 7r 2 + 3AC 2 fi 4 7 - 3\C 2 tt 4 + 2fi 4 7r 2 + 2m 4 vr 2 - 4ft 2 m 2 7r 2 )ln(^) ^ ^ 

128 7T 4 

m 2 

m 4 (46 2 7r 2 - 3AC 2 + 3AC 27 )ln(— ) 
A* 

128 " 7T 4 
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However, the nonlogariphmic term in Eq. (4.18), K(°\Q), includes a simple pole term: 

K(0) = (m 4 - (7 4 )[-8tt 2 (-1 + 7 )fe 2 - 3A(-1 + 7 ) 2 C 2 + 256£ 2 tt 4 ] + 47r 2 7 (m 2 - Q 2 ) 2 
^ ' ~ 256tt 4 

1 b 2 {Q 4 -m 4 ) 3 XC 2 {Q 4 -m 4 ) 



16 7T% 64 vr% 2 

(4.26) 

which can be made finite by the choice of D 2 as 

D = - J- A _ (4.27) 

2 16 7T% 64 7T% 2 ' 

Now the effective potential has been successfully regularized up to 5 2 order and is given 

by 

K ( ;}(0o) = ^-j (256A(7 2 4 7r 4 + 327r 2 (47r 2 (l + b 2 x) + 3xAC< 2 (ln(x/i) + 7 - 1))0 2 
+ 3A(7 2 (:r 2 ln(x/i) 2 - ln(/i) 2 ) + (6AC 2 (1 - 7 ) - 4vr 2 (2b 2 + l))HP) 

(4.28) 

+ ((67 - 6)Ax 2 C' 2 + 4vr 2 (l + 2b 2 x 2 ))\n(xft) + 3A( 7 - l) 2 (x 2 - 1)C 2 

+ 2(x - 1)7T 2 (46 2 ( 7 - l)(x + 1) + x - 3)) 

where we introduced dimensionless parameters x, p, , 0o defined by 

f2 2 = xm 2 , jj 2 = m 2 /fL, 0o = m 4>o , (4.29) 

and presented the effective potential in the units of m 4 . The effective potential includes two 
extra parameters Q and /i. In practical calculations the former may be determined by the 
principle of minimal sensitivity (PMS) [13] or alternatively, by the fastest apparent conver- 
gence (FAC) conditions ( see e.g. [11]) while the latter may be fixed by a renormalization 
scheme. 

Below we shall use the PMS which is natural for the renormalization scheme given in 
Eqs. (4.16). In accordance with the PMS the optimum value of Q will be given by the 
equation, dV e ^/dVl = , so, differentiating (4.28) with respect to x leads to the following 
'gap' equation: 

167r 2 x[46 2 7r 2 + 3AC' 2 ln(x) + 3AC 27 + 3AC 2 ln(/i)]0 2 

+ 3AC- 2 a; 2 (ln 2 (/i) + ln 2 (x)) + 27r 2 [2b 2 x 2 (2 7 - 1) + (x - l) 2 ] 

(4.30) 

+ [6XC 2 x 2 \n(jl) + (6x 2 A 7 - 3x 2 X)C 2 + 8x 2 6 2 7r 2 ]ln(x) 

+ x 2 (6AC , 27 - 3\C 2 + 86 2 7r 2 )ln(/i) + 3 7 x 2 A( 7 - l)C 2 = 



13 



which determines 11 as a function of 0o- Note that, this dependence, or, more exactly, 
dQ/d<f) should be explicitely used when calculating high order total derivatives of V^g with 
respect to <po, for example, in Eqs. (4.16). The first condition in Eqs. (4.16) gives: 

3A(7 2 (ln(//) - 1 + 7 ) + 46 2 tt 2 = (4.31) 

which determines hi as 

b = 3 AC^lnpZ) -I+7) (432) 

4 71-2 

and the second condition 

[2(7 + ln^)) 2 - 41n(/i) + 4 - 4 7 ]C 2 - 2( 7 + ln^)) 2 - 1 + 4 7 + 41n(/i) = (4.33) 
determines C2 as: 

, = 1 2( 7 + ln(/i)) 2 -41n( / i)-4 7 + l 

2 2 ( 7 + ln(/i)) 2 + 2-21n( / x)-2 7 ' 1 ' 1 

Now substituting this into (4.30) and considering the point 0o — 0, x — 1, (i.e. O = m at 
0o = 0) we get : 

2Aln 4 (/z) + 8A(-1 + 7 )ln 3 (/x) + A(-24 7 + 12 7 2 + ll)ln 2 (/i) 

(4.35) 

+ 2A(-1 + 7) (2 7 - 1)(2 7 - 3)ln(/x) + A(2 7 2 - 4 7 + 1)(-1 + 7 ) 2 = . 
One can easily check that this equation is satisfied for 

\n(fi) = \n(m 2 /n 2 ) = 1 - 7 (4.36) 
which leads to the following values for the finite counter terms: 

62 = , C 2 = ~. (4.37) 

Finally, the renormalized effective potential up to S 2 order is given by: 



(2) 1 [-3xAln(x) + 3xA+167r 2 ]0 2 
K//(0o) - "2^00 + 3^ 

87r 2 (a;- 1) -ln(x)(3x 2 A + 87r 2 ) 



(4.38) 



512vr 4 

with x being the solution of the gap equation, 

3xA[m(x) + 1] [xln(x) + 160 2 tt 2 ] - 4tt 2 (x - l) 2 = (4.39) 

Comparing Eqs. (4.38), (4.39) with the results of the "precarious" version of renormal- 
ization obtained by Stancu and Stevenson ( Eqs. (3.30) and (3.21) of [4]) , one can see 
that their mathematical forms are different. One of the reasons of this difference is that, 
in ref. [4] the finite part of two and three loop integrals ^3 and I4 were not given while in 
the present work we used the explicit form of the finite part ( see Appendix). The second 
reason may be hidden in the fact that, to make the effective potenial finite the most authors 
use the procedure "regularization after renormalization" while in the present work we first 
regularized the potential, i.e. made it finite, and then applied the renormalization conditions 
- Eqs. (4.16). Note that, the nessesity that the regularization should be carried out before 
imposing the gap equation was pointed out also in refs. [10,14]. 
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V. CONCLUSIONS 



We have developed a systematic method for computing the correction terms to the Gaus- 
sian effective action. They are given by the expansion of: 

exp{~[v 2 AW + v 3 A^ + v A A^]}exp {hjGj/2} (5.1) 

in powers of 5. Although our final result for unrenormalized PGEP is the same as that of 
refs. [3] and [4], the method proposed in the present work has some advantages. 

1) In refs. [3,4] the expansion parameter 5 is introduced directly into the Lagrangian at 
the very beginning of calculations. Instead, by using BFM we preliminary extracted the 
Gaussian part of the effective potential and then introduced an expansion parameter 5 in 
Eq. (3.17). This seems to be a more natural way, since in accordance with general rules of 
quantum mechanics exponent of any operator A should be understood as a Taylor expansion, 
i.e.: exp(i) = {1 + 5 A + (5A) 2 /2\ + (5Af/?>\ + . . .}\ s=1 . 

2) In the present approach Gaussian part of the effective action is extracted by introducing 
the primed derivatives (3.8). We have shown that, in our method the computation of 
correction terms to the GEP is greatly simplified. In fact, at each order of 5 it gives BFM 
diagrams (three in 5 2 order and six in 5 3 order) with a quadratic vertex v 2 = — 
f2g) + 6A O 2 , cubic vertex 4A O and quartic vertex A . These diagrams do not include a 
ring diagram (G xx I ) explicitly, which is included there through the relation SIq = Q 2 — 
12h\oIo(fl), giving rise to the ones (five in 5 2 order and 13 in <5 3 order) appearing in the 
methods used in ref.s [3,4]. 

3) In the present method it is easier to solve equations (2.7) -(2.9) of BFM than similar 
equations in the method used in [4]. In fact, in the method of refs. [3,4] one has to solve the 
equations T[0] = W[J] — J(f), <fi = 5W/5J to define and J at each order of 5. Although the 
corresponding solutions are simple in 5 2 order 4 , J = 5(vi + 3t>o/o) + 0(S 2 ), they are very 
complicated at higher orders. On the other hand, we have shown that due to the advantages 
of BFM one may simply put J = at each order of 5. Thus using the method proposed in 
the present paper it will be much simpler to calculate higher order corrections to GEP in 
general. 

We have shown that the variational perturbation expansion of 4 scalar theory in four 
dimensions can be successfully renormalized by introducing appropriate counter terms in 
such a way that the whole Lagrangian has the same polynomial form as the bare one. 
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4 Eq. (2.29) of ref. [4]. 
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Appendix 



APPENDIX: EXPLICITE EXPRESSION FOR DIVERGENT INTEGRALS. 

Here we bring explicit expressions for the divergent integrals defined as : 

d A p If d 4 p\n{p 2 + fl 2 ) 



J0(fi) -/W(^' m) ~2S (27T) 

J 2 (fi) _ f d*pG 2 (p) 2 2 

— = J (2vr)4 ' G(p) = l/(p +0 ), 



I 3 (Q) _ f d*pd*q G(p)G(q)G(p + q) 



A 



(Al) 



/ 



3! ~J (2tt) 8 
J 4 (fi) _ /■ d*pd A qd A k G{p)G(q)G(k)G(p + q + k) 



-/ 



4! 7 (2tt) 



12 



in four dimension. Integrals Jo(fi),li(fi) and l2(ft) are quite simple and may be found elsewhere 
[8]. 



, ^ 2 [-l + 7o + ln(^^)] 



8vr 2 e 16vr 2 
1 ft 2 

^ = 4^ " 8 ^ ( A2 ) 

h(n) = h{m) + -(Q 2 - m 2 )h{m) - -(ft 2 - m 2 ) 2 I 2 {m) 
2 8 

^ m 4 [2x 2 ln( x) - 2x + 2 - 3(x - l) 2 ] 



128vr 2 

where x = £l 2 /m 2 , e is defined as d = 4 — e , d- space - time dimension, 70 = 0.5777215 and \x 
is an arbitrary constant with mass dimension which usually appears in dimensional regularization 
scheme. 

The two loop and three loop integrals h(Cl) have been accurately calculated in refs. [15] 

including finite parts: 
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i n 2 

3 ,L1 -2470 - 241n (i -,) + 36 



= -^-l^ 2 24-, + + 121n 2 ( 



2567T 4 "" \~V e v 4vr^ 2 " 



2 



+ (24 7 o-36)ln(i-^)+ 1.2 A - 3(> 7l + L2~, 2 - ;r • • 12 



92 , „ 4 ft 2 
- 247o - 241n(- 



hm = ^ 1 _3 4tt/^ 

4V ; 512vr 6 e 3 e 2 



1 Q2 i q2 

35 - 46 7o + 18 7 2 + tt 2 + (36 7o - 46)ln(-— ) + 181n 2 (-- 



k 4 7T/X 2 4 7T / [/ 2 ' 



, ,1 SI 2 w 105 „ , 3 9N , ? .l n 2 . .69 . ,0,1 ft 2 

+ In(-— .)(-- - + 69 7o - 27 7o 2 - -tt 2 ) + ln 2 (-^)(- - 27 7o ) - 91n 3 (- 



v 4 7r/i 2 ,u ,u 2 7 v 4tt/x 2/v 2 ,uy Mvr^ 2 

where A = —1.171953. Note that in the main part of the text we used notation 
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FIGURE CAPTIONS 



Fig. 1(a). BFM Feynman diagrams contributing to the effective potential in 5 2 order (see 
Eq. (3.23)). Solid lines represent the propagator G(p) = l/(p 2 + fi 2 ). The vertices 
v 2 = (m^ — Qq)/2 + 6A O 2 , v 3 = 4A O and t> 4 = A are marked by diamonds, circles 
and squares, respectively. 

Fig. 1(b). Feynman diagrams contributing to the effective potential in 5 2 order in the 
method of refs. [3,4] (see Eq. (4.1)). These are obtained by an attachment of the 
ring diagram, Iq , denoted here by the small circle, to the each two body vertex v 2 of 
the FIG. 1(a). Solid lines represent the propagator G(p) = l/(p 2 + fi 2 ). The vertices 
U2 = (ml — Qq)/2 + 6Ao0 2 , t>3 = 4Ao0, t>4 = Ao and w 4 = 6Ao are marked by diamonds, 
circles, filled and opened squares, respectively. 

Fig. 2(a). The same as in FIG. 1(a). but in 5 3 order. 

Fig. 2(b). Feynman diagrams contributing to the effective potential in 5 3 order of ref. [3]. 
These are obtained by attachment of a ring diagram, Jo , denoted here by the small 
circle, to the each two body vertex f 2 of the Fig. 2 (a). The notations are the same as 
in FIG. 1(b). 

Fig. 3. The same as in FIG. 1(a) but in 5 A order. Note that, only those, including quadratic 
vertex i> 2 are presented. Implementation of a ring diagram to each t> 2 vertex with all 
possible ways will give 18 additional diagrams. 
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FIG. la 
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FIG. lb. 
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FIG.2a. 
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FIG.3. 
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